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Abstract 

The analysis of the solution with full replica symmetry breaking in the vicinity of 
Tc of a general spin glass model with reflection symmetry is performed. The leading 
term in the order parameter function expansion is obtained. Parisi equation for the 
model is written. 
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The crucial role of the reflection symmetry for the character of phase transition 
in nonrandom mean-field (MF) models is well known (see, e.g., the textbook 
[1]). Generally speaking the presence of the terms without reflection symmetry 
causes the first order phase transition, while in the absence of such terms the 
transition is of the second order. Usually this result is obtained in the frame of 
the phenomenological approach based on the Ginzburg-Landau (GL) effective 
Hamiltonian that can be easily obtained for any Hamiltonian through the 
Hubbard-Stratanovich identity for the partition function. 

In the case of random MF models the behavior of systems with reflection 
symmetry is quite different from those without it. In Ref. [2] this problem 
was investigated in the frame of replica symmetric (RS) approach. The role 
of the reflection symmetry for the behavior of RS solution for spin-glass-like 
MF systems was considered and a kind of symmetry rule for the type of the 
growing of glass order was formulated: if in the nonrandom (pure) system the 
transition to the ordered phase is of the second order, then in the corresponding 
random system the glass regime appears as a result of a phase transition; if 
the transition in the pure system is of the first order, then in the random 
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system the glass order parameter grows continuously on cooling. In fact, one 
can imagine that both the first order phase transition in nonrandom systems as 
well as the continuous growing of the glass order in random systems are caused 
by some kind of internal fields appearing due to the algebra of operators U. In 
this case we have a close analogy with Sherrington-Kirkpatrick (SK) model 
in an external field. 

The scenarios of replica symmetry breaking (RSB) for the mentioned two cases 
differ essentially one from another because the absence of refiection symmetry 
results in a special form of RSB free energy functional for random case, too. 
So, the absence of reflection symmetry can lead to continuous growing of the 
Edwards-Anderson order parameter as well as to a jump at the transition to 
one stage replica symmetry breaking (IRSB) phase and other specific features 
(see, for example, Refs. [3,4,5,6]). In the refiection symmetrical random sys- 
tems the glass regime appears as a result of a bifurcation from the trivial solu- 
tion for the order-parameter function. The prototype model is the well known 
Sherrington-Kirkpatrick model [7], the fully connected Ising spin model with 
quenched random exchange magnetic interactions. This model is described by 
a full replica symmetry breaking (FRSB) Parisi solution [8,9,6,10]. 

Although the features of SK model have been thoroughly investigated , com- 
paratively little attention has been paid to the analysis of natural generaliza- 
tions of SK model (for example the systems of higher spins). So, we believe 
that it would be quite useful to consider a generalized spin glass model where 
Ising spins are replaced by arbitrary diagonal operators JJ with TtU = 0. In RS 
approximation this model was considered in [2] . In the present paper we con- 
tinue the investigation of this model and consider in detail replica symmetry 
breaking in the case when the system possesses refiection symmetry. 

Let us consider a system of particles on lattice sites i, j with Hamiltonian 

h^~Y.JAUj, (1) 

where U is an arbitrary diagonal operator with Trll — 0, Jy are quenched 
Gaussian interactions with zero mean 

with J = J/y/N. 

Using replica approach (see, e.g. [7]) we can write the free energy averaged 
over disorder in the form [2] 
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Here t = J/kT. The saddle point conditions for the free energy give the glass 
order parameter 



and auxiliary order parameter 
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As one index quantity are usually replica invariant (see e.g. [5] and [6]) we 
put — p from the very beginning. 

In this paper we study the model (1) for operators U with zero trace for all 
odd powers: 

Tr [;7(''=+^)] = (7) 

for all integer k. In the replica symmetric (RS) approximation the free energy 
(3) has the form: 



F = -NkTit'i-t^^+ 
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Here 



The extremum conditions for the free energy (8) give the following equations 
for the glass order parameter: 



and an auxiliary equation 
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Here 



^ dz I 



If the condition (7) is fulfilled then Eq. (9) has a trivial solution for the glass 

order parameter q^^s = at any temperature. The nontrivial solution appears 
at the point Tc defined by the bifurcation condition for the Eq. (9): 



1 - tlp{tc) = 0. 



(11) 



The replica symmetric solution is stable unless the replicon mode energy A is 
nonzero [11]. For our model (1)- (3) we have: 
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In our case (7) the condition Xrepi(RS) = is equal to the equation (11) 
due to the fact that it is the trivial solution that bifurcates. Performing one 
stage rephca symmety breaking (IRSB), then 2RSB, 3RSB etc. we can see 
that the r.h.s. of the equations for the glass order parameters always contains 
TriUexp{e„RSB)] _ Therefore at each stage of RSB one possible solution is a trivial 
one and we can consider the appearance of nRSB solution as the bifurcation 
from the trivial (nRSB) solution. In such the equation X(nRSB) = 

coincides with the corresponding bifurcation condition (11). For example, in 
the case of 2RSB we have 

Kepl{lRSB) = 1 — 

i'/ fr^i^^Ja ^77^ / "^y"" (Tr[exp^iRSB(-2,2/)])' 
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Here 
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The IRSB solution is characterized by m and two vahics r and r + v of q°'^. 
For r — and v — the equation (13) reduces to Eq. (11). 

This means that the RSB solutions of different stages at least may exist at 
temperature T < T^, defined by the condition (11) and encourage us to look 
for the full replica symmetry breaking (FRSB) [8,9]. 

As an example let us consider the Hamiltonian (1) with U = Q + r}V, Q ~ 



3Jz — J{J + 1),V — V^iJx ~ Jy )■ The model describes random quadrupole 
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Fig. 1. IRSB order parameters and Ai.epi(iRSB) for the case i] = l/\/3 (analogue for 
spin-one glass model). Curves for p^s and Pirsb coincide. 



interaction. This model has been considered for J = 1, = 1, 0, — 1 in Ref. [2] 
and for J = 2 in [12] in RS approximation. Here we restrict ourselves to the 
case J = 1. It easy to show that = 2 - Q, = 2 + Q, QV = VQ = V 
and 

(g + vVf = 2(1 + r?2) + 2riV - (1 - ri^)Q. (14) 
(Q + vVf = 2{?>rf - 1) + 3(1 + r/2)(Q + r,V)]- 

i:i{Q + r]Vf = 6(377^ - 1). 

Here > is a tuning parameter. The behavior of the system is quite different 
in the cases r] = l/y/S (analog of spin-one glass model [15]) and r] ^ l/^/S. 

In figures we present the behavior of the IRSB order parameters and Xrepi{iRSB) 
for this case (77 = 1/ -\/3) and for the spin-two glass model. 

In general case in FRSB scheme the series expansion in small 6q°'^ (up to 
fourth order) of the free energy (3) in the vicinity of has the following 
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Fig. 2. IRSB order parameters, Ai.epi(RS) and Ai.epi(iRSB) for spin-2 glass model. 
Curves for p^s and Pirsb coincide. 
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The indices in (15) are different only if they are connected through Sq"^^ . We 
beheve that the parameter p (5) changes slowly and is equal to Pc in the 
vicinity of tc- 

Now using Parisi rules [8,9] we can write the free energy (15) in the continuous 
case of FRSB in the form: 
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The usual approach to find the extremum of the free energy with variation 
with respect to q{x) and further differentiations with respect to x leads in the 
main approximation to 



q{x) = < 



cx, < X < ^; 
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c — 



(18) 



In close analogy to the procedure described in [13,14] one can easily obtain 
the expansion of the free energy and of the order parameter q{x) up to higher 
orders. One can use the approach of Parisi equation [8] for obtaining numerical 
solution. In our model the free energy in an external field h can be written in 
the form 
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NkT 4 
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where 0(0, y) is the solution at x = of the Parisi equation 



q[x) 
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(20) 



with the boundary condition 

0(l,y) = ln'iv|exp 



tyU+\^{p-q{l))U^ 



(21) 



The derivative with respect to x are denoted by overdots and the derivatives 
with respect to y by primes. In obtaining (20) - (21) the fact that Tr = 
was not used. However we assume that the FRSB solution exists and g(0) = 
0. 



The results (16)- (21) reduce to well known form in the case of Sherrington- 
Kirkpatrick model [7] when p = = 1. As another example the case of 
spin-one can be considered. The replica symmetric solution for this model was 
given in [15]. In this case "ip — 2-\-Q where Q is the quadrupole moment. The 
order-parameter function in FRSB scheme was obtained for spin-one model 
in [16] (i^ = 0). Our results (17) - (18) can be reduced to those of [16]. 

Now we should like to emphasize the fact that if the Hamiltonian of a model 
is the sum of different Hamiltonians, each of them possessing reflection sym- 
metry, the behavior of the system can be completely different relative to that 
discussed here in the case of nonzero mixed terms of third order. This can be 
demonstrated using the toy-model with the Hamiltonian 



(22) 



where are quenched Gaussian interactions with zero mean and 



C/(i) = 
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Tr(f/('))3 = 0, but Tr [(?/(') )2f/(^)] ^ if / ^ A; and this results in the appearance 
of the term Y^{5q'^^Y the free energy expansion [17]. The model described 
by the Hamiltonian (22) is isomorphic to three state Potts model and has an 
IRSB solution which is stable in a certain temperature interval as was shown 
in [18]. 

To conclude, in the present paper we have shown the possibility of Parisi full 
replica symmetry breaking for rather general model with reflection symmetry. 
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The model is a generalization of SK model. We obtain the main term of the 
expansion of the FRSB solution for the spin glass order parameter near T^. 
The Parisi equation is also written in the general case. 
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